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I. INTRODUCTION 

One of the remaining problems in low energy QCD 
is the quantitative field theoretical description of the 
confinement-deconfinement phase transition. Apart from 
its genuine importance for a first principle understanding 
of the confining physics in QCD it also is a key input for 
the evaluation of the QCD phase diagram. 

In the past decade much progress has been achieved 
both in continuum studies as well as with lattice compu- 
tations for our understanding of the low energy sector of 
QCD, for reviews see e.g. [iTIl, S 0, Hi- For an analyt- 
ical description of the low energy sector, topological de- 
grees of freedom are likely to play an important role for 
the confinement-deconfinement phase transition as well 
as for chiral symmetry breaking, see e.g. Q. The lat- 
ter has been very successfully described within instanton 
models, whereas the confining properties of the theory 
are harder to incorporate within semi-classical descrip- 
tions. Indeed, tracking down those topological degrees of 
freedoms relevant for confinement in the physical vacuum 
has its intricacies as the physical vacuum is more likely to 
contain a rather dense packing of topological configura- 
tions, making their detection difficult. Moreover, models 
of confinement are rather based on topological defects 
instead of stable topological objects, the construction of 
which out of these defects is plagued by non- localities. 

Still, these defects are manifest in the Polyakov loop, 
the order parameter in pure Yang-Mills theory 0], and 
can be extracted by an appropriate gauge fixing, see e.g. 

Q . Gauge fixing is also mandatory in most contin- 
uum formulations of QCD for removing the redundant 
gauge degrees of freedom. This is mostly seen as a lia- 
bility of such an approach, as a formulation of QCD in 
gauge-variant variables complicates the access to gauge 
invariant observables. However, gauge fixing is nothing 
but a reparameterisation of the path integral and can be 
used for even facilitating the computation of at least a 
subset of observables. Indeed, this point of view has been 
exploited much in the discussion of confinement mecha- 
nisms based on topological defects. More recently it also 
has become clear that these are not competing physics 
mechanisms but rather different facets of the same global 



physics picture which still awaits a fully gauge invariant 
description, see e.g. [13] . Despite this final step we have 
learned much from the combined investigations which to- 
gether built a nearly complete mosaic. 

The effective potential of the Polyakov loop has also 
been used as an input for effective field theories that 
give some access to the QCD phase diagram At 
finite temperature and vanishing density, these models 
have led to impressive results in particular for thermo- 
dynamical quantities. At finite chemical potential, the 
back-reaction of the matter sector to the gauge sector is 
difficult to quantify in these models, and the chiral and 
confinement-deconfinement phase transitions are sensi- 
tive to the details of this back-reaction. This also holds 
true for the question of a quarkyonic phase with confine- 
ment and chiral symmetry at finite density p^ . For an 
extension of these models one has to resort to a field- 
theoretical description of the gauge sector which allows 
to systematically study the impact of a finite chemical 
potential on the dynamics of the gauge field, [l^ . 

In summary, the evaluation of Green functions of the 
Polyakov loop allows for a direct access to the physics 
in the strongly coupled sector of QCD, and in particu- 
lar the confinement-deconfinement phase transition. In 
the present work we initiate a non-perturbative study 
of QCD in Polyakov gauge. In this gauge the Polyakov 
loop takes a particularly simple form and is directly re- 
lated to the temporal component of the gauge field. After 
integrating-out the spatial components of the gauge field, 
and formulated with Polyakov loop variables, the gauge 
field sector of QCD resembles a scalar model. The dy- 
namics of low energy Yang-Mills theory is then incorpo- 
rated by evaluating Wilsonian flows for the effective ac- 
tion [3, Q, [ij, U^, [3, II3l ■ We derive the flow equation for 
QCD in Polyakov gauge, and solve it for the full effective 
potential of the Polyakov loop. Due to the formulation in 
Polyakov gauge a simple truncation already suffices to en- 
code the physics of the confinement-deconfinement phase 
transition. The results include the temperature depen- 
dence of the Polyakov loop, and the critical temperature. 
We also compare the present approach to lattice studies 
p^ . and to a recent continuum computation in Landau 
gauge dil. 
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II. QCD IN POLYAKOV GAUGE 

In QCD with static quarks the expectation value of 
a static quark serves as an order parameter for 

confinement. It is proportional to the free energy Fq of 
such a state, {q{x)) ~ exp(— /3Fq), where /3 = 1/T is the 
inverse temperature. Hence in the confining phase at low 
temperature, the expectation value vanishes, whereas at 
high temperatures in the deconfined phase, it is non-zero. 
The Polyakov loop variable, 0, is the creation operator 
for a static quark, 

L(f) = _LtrP(f), (1) 

where the trace in ([T]) is done in the fundamental repre- 
sentation, and the Polyakov loop operator is a Wegner- 
Wilson loop in time direction, 

'P(x) = "P exp ^i^y dxoAoixo,x)^ ■ (2) 

Here V stands for path ordering. We conclude that 
{q{x)) ~ {L{x)), and thus the negative logarithm of the 
Polyakov loop expectation value relates to the free en- 
ergy of a static fundamental color source. Moreover, (L) 
measures whether center symmetry is realised by the en- 
semble under consideration, see e.g. [ll. [gLItIIsI. [qI. [Toj . 

More specifically we consider gauge transformations 
U{xo,x) with U{0,x)U~^{l3,x) = Z, where Z is a center 
element. In SU{2) the center is Z2, whereas in physical 
QCD with SU{3) it is Z3. Under such center transforma- 
tions the Polyakov loop operator 7^(x) in ^ is multiplied 
with a center element Z, 

V{x)^ZV{x), (3) 

and so does the Polyakov loop, L{x) Z L{x). Hence, 
a center-symmetric confining (disordered) ground state 
ensures (L) = 0, whereas deconfinement (L) 7^ signals 
the ordered phase and center-symmetry breaking, 

T<T,: {L{x)) = , = 00 , 

T>r, : (L(f))^0, <oo. (4) 

The expectation value of the Polyakov loop can be de- 
duced from the equations of motion of its effective po- 
tential Vl[{L)]. We shall argue, that the computation of 
the latter greatly simplifies within an appropriate choice 
of gauge. Indeed, gauge fixing is nothing but the choice 
of a specific parameterisation of the path integral, and 
a conveniently chosen parameterisation can simplify the 
task of computing physical observables. 

In the present case our choice of gauge is guided by 
the demand of a particularly simple representation of 
the Polyakov loop variable (H]). A gauge ensuring time- 
independent ^0 leads to both, a trivial integration in ^ 
and renders the path ordering irrelevant. Having done 



this one can still rotate the Polyakov loop operator P{x), 
(HI), into the Cartan subgroup. The above properties are 
achieved for time-independent gauge field configurations 
in the Cartan subalgebra, i.e. Ao{to,x) = Aq{x). For 
SU{2) this reads 

Aoito,x)=Aoix)^ (5) 

and entails a particularly simple relation between ^0 and 
L, 

L{x) ^ cos ]^gf3Ao{x) , (6) 

Note that this simple relation is not valid on the level 
of expectation values of L and Aq, in SU{2) we have 
(L) 7^ cos ^gf]{AQ). However, in the present work we 
consider an approach that gives direct access to the ef- 
fective potential V;,ff[(j4o)] for the gauge field, as distin- 
guished to those for the Polyakov loop, Ucff[{L)] ^. 

Here, we argue that L[(j4o)] also serves as an order 
parameter: to that end we show that the order parame- 
ter (i[Ao]) is bounded from above by L[{Ao)]. It follows 
that L[(Ao)] is non- vanishing in the center-broken phase. 
Furthermore we show that in the center-symmetric phase 
with vanishing order parameter, {L[Ao]) = 0, also the ob- 
servable L[(74o)] vanishes. For the sake of simplicity we 
restrict the explicit argument to SU{2), but it straight- 
forwardly extends to general SU{N). First we note that 
we can use ([SJ for expressing the expectation value of Aq 
in terms of L, 

^.9/3(^0) = (arccosL) . (7) 

We emphasise that the rhs of ^ defines an observable as 
it is the expectation value of an gauge invariant object. 
This observable happens to agree with (Aq) in Polyakov 
gauge. It follows from the Jensen inequality that the 
expectation value of the Polyakov loop, the order param- 
eter for confinement, is bounded from above by L[{Aq)], 
see 

L[{A„)] > {L[Ao]) . (8) 

for gauge fields gf3{Ao)/2 e [0, 7r/2]. Note that it is suf- 
ficient to consider the above interval due to periodicity 
and center symmetry of the potential. This means we 
restrict the Polyakov loop expectation value to (i) > 0. 
Negative values for (L) are then obtained by center trans- 
formations, L ±L. Eq. ([5|) is easily proven for SU (2) 
from ([6]) as cos(a;) is concave for x £ [0, 7r/2]. Thus, for 
(L) > it necessarily also follows that g(3{AQ)/2 < 7r/2. 

In turn we can show that g/3{Ao)/2 = 7r/2, if the 
Polyakov loop variable {L[Aq]) vanishes. This then en- 
tails that L[(Ao)] = 0. To that end we expand L about 



^ A reformulation in terms of the Polyakov loop variable only along 
the lines outlined in [TtII is in progress. 
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its mean value (L), that is L = (L) + 5L. Inserting this 
expansion into ([7]) we arrive at 



^gf3{An) = arccos(i) - ^ 



= {6L)+0{{6L')) 



(9) 

In the center-symmetric phase (L) = 0, c.f. ([4]). Un- 
der center transformations L transforms according to ([3|) 
L Z L with Z = ±1 and hence SL Z SL. It follows 
that ((5i2"+i) = Z((5L2«+1) ^ 0, and all odd powers in 
^ vanish. The even powers vanish since arccos is an 
odd function and hence has vanishing even Taylor coeffi- 
cients arccos*^^"-* (0). Thus, in the center-symmetric phase 
we have 



1 TT 

-gP{Ao) = arccos(L) = - . 



(10) 



In summary we have shown 

T<T,: L[(Ao)]^0 ^ ig/3(Ao(x)) = | , 

T>T,: L[{Ao)]^0 ^ ig/3(Ao(f)) < | .(11) 

We conclude that (Aq) in Polyakov gauge serves as an or- 
der parameter for the confinement-deconfinement (order- 
disorder) phase transition, as does L[{Ao)]. Thus, we 
only have to compute the effective potential T4ff[(^o)] in 
order to extract the critical temperature, and e.g. critical 
exponents. This potential is more easily accessed than 
that for the Polyakov loop. It is here were the specific 
gauge comes to our aid as it allows the direct physical in- 
terpretation of a component of the gauge field. This prop- 
erty has been already exploited in the literature, where 
it has been shown that (Aq) in Polyakov gauge is sen- 
sitive to topological defects related to the confinement 
mechanism |9|. 



III. QUANTISATION 

We proceed by discussing the generating functional of 
Polyakov gauge Yang-Mills theory. For its derivation we 
use the Faddeev-Popov method. Specifying to SU{2), the 
Polyakov gauge ([5]) is implemented by the gauge fixing 
conditions 

9otr asAo = , tr (ai ± ia2)Ao ^ , (12) 

where the ct^ are the Pauli matrices. However, the gauge 
fixing is not complete. It is unchanged under time- 
independent gauge transformations in the Cartan sub- 
group. These remaining gauge degrees of freedom are 
completely fixed by the following conditions, 

di J dxoti a^Ai =0, 82 J dxodxitr a3A2 = 0, 

^3 / dxodxidx2 tr (T3A3 — . (13) 



The gauge fixings ([T3|) are integral conditions and are 
the weaker the more integrals are involved. Basically 
they eliminate the corresponding zero modes. This can 
be seen directly upon putting the theory into a box with 
periodic boundary conditions, T^, see e.g. @. 

The gauge fixing conditions (fT^ . (|13p lead to the 
Faddeev-Popov determinant 



Afp[A] - (2T)2 



V 2T 



(14) 



the computation of which is detailed in appendix 1X1 The 
integration over the longitudinal gauge fields precisely 
cancels the Faddeev-Popov determinant in the static ap- 
proximation diAg = 0, see Appendix Finally we are 
left with the action 

SMA] ^ -1/3 [ d'xAod^Ao (15) 



\ I d^xA'', 



2\ab |_ q2 ^ab 



with D'^ = d^i""^ + A^gJ""^^ and transversal spatial gauge 
fields, diA^^i = 0. The generating functional of Yang- 
Mills theory in Polyakov gauge then reads 

Z{J\ = J dAadA^,, cxp{-S,ff[A] 

+ j d^xJoAo + J d^xJ^,,A^,,j . (16) 

In (jl6p we have normalised the temporal component Jq 
of the current with a factor /3. The classical action Ses 
is inherently non-local as is contains one-loop terms, the 
Faddeev-Popov determinant as well as the integration 
over the longitudinal gauge fields. 

Instead of computing Z[J] in (fT6)) we shall compute 
the effective action T within a functional rcnormalisation 
group approach [1, i, [ll, [H, [H, [l3. To that end we 
introduce an infrared cut-off for the transversal spatial 
gauge fields and in the temporal gauge fields by modify- 
ing the action, S Ses + ASkiAo] + AS±,k[A±], with 
infrared scale fc, and cut-off terms 



ASk[Ao] 



1 



(3 / d^xAoRo,kAo 



AS, 



k,±[ 



d^xA'i.R^^kA'i, 



(17) 



The regulators Rk in ^17} arc chosen to be momentum- 
dependent and required to provide masses at low mo- 
menta and to vanish at large momenta. For fc ^ they 
vanish identically. 

They can be written as one single regulator Ra.^v, 
which is a block-diagonal matrix in field space with en- 
tries i?y!i,oo = RQ,k and RA,ij = R±.k^±,ij, where the 
transversal projector is defined by 



n 



(18) 
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The above structure is induced by the fact the A±^i are 
transversal, and hence Rj_^k only couples to the transver- 
sal degrees of freedom. 

The flow of the cut-off dependent effective action Tk is 
then given by Wetterich's equation [3, [I^, [1^1 for Yang- 
Mills theory d, in Polyakov gauge, 



momentum dependence 



2 J (2^)3 \t^)+R^ 



T y-^ f (Pp 



(27r)3 \rf) 



dtR±,k , (19) 



where t is the RG time t = \n{k/A), and A is some refer- 
ence scale. 



IV. APPROXIMATION SCHEME 



Eq. p9|) together with an initial effective action at 
some initial ultraviolet scale k = Ayv provides a defini- 
tion of the full effective action at vanishing cut-off scale 
A: = via the integrated flow. For the solution of ([TO]) 
we have to resort to approximations to the full effective 
action. In gauge theories such an approximation also re- 
quires the control of gauge invariance, see e.g. (l7j . 

Here we shall argue that in Polyakov gauge a rather 
simple approximation to the full effective action already 
suffices to describe the confinement-deconfinement phase 
transition, and, in particular, to estimate the critical tem- 
perature. We compute the flow of the effective action 
r[^Oi^±] in the following truncation 



lim (L(O)Lt(p)), 



1 



Pa 



lino ^2 ((^^)2+p2)2 



1 1 



7r2 (/3cr)3 ■ 

(22) 

We conclude that the Polyakov loop variable has a mas- 
sive propagator. This directly relates to a massive prop- 
agator of ^0 in Polyakov gauge. 

The approximation scheme is fully set by specifying 
the regulators i?o,fc and R±^k- Naively one would iden- 
tify the cut-off parameter k in the regulators with the 
physical cut-off scale fcphys- For general regulators this 
is not possible and one deals with two distinct physi- 
cal cut-off scales, /co^phys and A:^^phys related to i?o,fc and 
R±^k respectively, for a detailed discussion see How- 
ever, within the approximation (|20p it is crucial to have 
a unique effective cut-off scale fcphys = fco,phys = fc_L,phys, 
as different physical cut-off scales A;o,phys 7^ k±, 

phys necGS 

sarily introduce a momentum transfer into the flow which 
carries part of the physics. This momentum transfer is 
only fully captured with a non-local approximation to 
the effective action. In turn, a local approximation such 
as (PO)) requires fco.phys = ^"i.phys- In other word, a local 
approximation works best if the momentum transfer in 
the flow is minimised. More details about such a scale 
matching and its connection to optimisation (TtI . |20| can 
be found in • Note in this context that in the present 
case wc also have to deal with the subtlety that Aq only 
depends on spatial coordinates whereas Aj_ is space-time 
dependent. However, the requirement of minimal mo- 
mentum transfer in the flow is a simple criterion which 
is technically accessible. 

More specifically we restrict ourselves to regulators (2]| 



rk[Ao,Aj 



d^3 



Zn 



Aod^Ao + Vk[Ao] 



^J^d^xZ,Al 



+ Al, (20) 



with fc-dependent wave function renormalisations Zq, Zi. 
The effective action (PU)) relates to the order parameter 
{L{x)) as well as its two point correlation {L{x)L'^ (y)) 
via the effective potential Vcfti^o] = Vfc[Ao] as explained 
in sectionini The expectation value {L{x)), or L[(Ao)], is 
used to determine the phase transition temperature as 
well as critical exponents. The temperature-dependence 
of the Polyakov loop two-point function relates to the 
string tension. In the confining phase, for T < Tc, and 
large separations |af — ?/| oo, the two-point function 
falls off like 

lim (L(f)Lt(y)),~exp{-/3a|x-y|} . (21) 

a;— y I — >C50 

Here, {■ ■ ■)c stands for the connected part of the related 
correlation function, i.e. {L{x)L'^ [y]) c = {L{x)L^ (y)) — 
{L{x)){L(y)). In turn, its Fourier transform shows the 



^A,00 — ZoRopt,k{p ) , RA,ij — ZiU±^ij{p)Ropt,k±{'P ) , 

_ (23) 
where [201 



RoptAf) = ik-p'nk'-f) 



(24) 



The detailed scale-matching argument is deferred to Ap- 
pendix [Cl and results in a relation k± = k±{k) depicted 
in Fig. [7] in the appendix. It is left to determine the 
effective cut-off scale fcphys- This cut-off scale can be 
determined from the numerical comparison of the flows 
of appropriate observables: one computes the flow with 
the three-dimensional regulator Ropt.k^ip^) in as 
well as with the four-dimensional regulator i?opt. fcphys (p^)- 
Then the respective physical scales are identified, i.e. 
k±.phys{k±_) = fcphys- The results for this matching 
procedure are depicted in Fig. [5] in Appendix [Cl An- 
other estimate comes from the flow related to the three- 
dimensional v4o-fluctuations, where we can directly iden- 
tity fcphys = fc- We use the above choices as limiting cases 
for an estimate of the systematic error in our computa- 
tion. Our explicit results are obtained for the best choice 
that works in all physics constraints. 
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V. FLOW 

We are now in the position to integrate the flow equa- 
tion (Uni). To begin with, we can immediately integrate 
out the spatial gauge fields A± for Zi = 1, that is the 
second line in (fTQ]) . This part of the flow only carries 
an explicit dependence on the cut-off k, details of the 
calculation can be found in Appendix [B] It results in 
a non-trivial effective potential Vl^fc[Ao] that approaches 
the Weiss potential [23 in the limit k/T — > 0, and falls off 
like exp(— /3fcx(fc)) cos{gPAo) in the UV limit k/T 00, 
see Fig. [T] In terms of the effective action, after the in- 
tegration over A±, we are led to an effective action of 

TklAo] =p j <fx {^^{dA^f + A14[Ao] + V^^k[Ao_ 

(25) 

Eq. (123 follows from ^ with Tk[AQ\ = Vk[Ao,A±_ = 0], 
and 



Vu[Ao]^^VM + V^^k[Aa] 



(26) 



The full effective potential is given by V^ff[Ao] = 
AVfc=o[^o] + ^±,/c=o[^o]- Wc are left with the task to 
determine AV^, which is the part of the effective po- 
tential induced by Ao-fluctuations. In Polyakov gauge, 
these fluctuations carry the confining properties of the 
Polyakov loop variable, whereas the spatial fluctuations 
generate a deconfining effective potential for Aq, see Ap- 
pendix[Bj We emphasise that this structure is not present 
for spatial confinement which is necessarily also driven 
by the spatial fluctuations, and solely depends on these 
fluctuations in the high temperature limit. We hope to 
report on this matter in the near future. 

Here we proceed with the integration of the flow for 
the potential AT4. To that end we reformulate the flow 
(fT9|l as a flow for AVfe with the external input Vj^^fc, see 
The flow equation for AVfe reads 



dtRo.k 



(2^)3 Zop' 



dX{AVk + Vj_,k) + Ro.k 



(27) 

With the specific regulator Rk in we can perform the 
momentum integration analytically. We also introduce 
the scalar field ip = gPAo, and arrive at 



/39fcA14 



(l + 77o/5)fc2 



3(27r) 



(28) 



where the coupling gl has to run with the effective cut-off 
scale A:pijys, and is estimated by an appropriate choice of 
the running coupling a^, 



2 _ 9 



ith 



47ras(fc^ 



phys 



(29) 



see also Appendix [Bj This also entails that the anoma- 
lous dimension rjo is linked to the running coupling by 



-5tloga^(fcphys) ■ 



(30) 



At its root (|28|) is an equation for the dimensionless ef- 
fective potential V = P'^Vk in terms of V± = P'^V±^k and 
AV = P^AVk- The infrared RG-scale k naturally turns 
into the modified RG-scale k = k(3, that is all scales are 
measured in temperature units. Then the fiow equation 
is of the form 



dr.AV 



(l + 77o/5)fc2 



3(2^)2 



l + ^dUV±+AV) 



(31) 



The potential V and hence AV has a field-independent 
contribution which is related to the pressure. For the 
present purpose it is irrelevant and we can conveniently 
normalise the fiow PT|) such that it vanishes at fields 
where d^{V± + AV) = 0. This is achieved by subtracting 
2(1 + rio/5) fc2/(3(27r)2) in ^ and we are left with 



9^Al/=- — 



5 



1 



gl dl (y, + AV) 
AV) 



f.dl (Fx 



(32) 



where wc have kept the notation df,AV for df.AV ~2{ \ + 
770/5) k^ / (3(27r)2). In this form it is evident, that the flow 
vanishes for fields where ^'^{y^_-\- AV) = 0, i.e. once a re- 
gion of the potential becomes convex, this part is frozen, 
unless the external input V± triggers the flow again. 

We close this section with a discussion of the quali- 
tative features of It resembles the flow equation 
of a real scalar field theory, and due to Vj_, the flow is 
initialised in the broken phase. It relies on two external 
inputs, Vx and a^. 

The first input, Vj_, is computed in a pcrturbativc ap- 
proximation to the spatial gluon sector, and its computa- 
tion is deferred to Appendix [Bl It is shown in Fig. [1] for 
various values of the RG time fc, and approaches the pcr- 
turbativc Weiss potential [2^ for vanishing cutoff fc = 0. 
Wc have argued that within Polyakov gauge this approx- 



<> 




FIG. 1: Vj_ for different values of k 

imation should capture the qualitative feature of its con- 
tribution to the Polyakov loop potential. We emphasise 



6 



again that this is not so for the question of spatial con- 
finement, and the related potential of the spatial Wilson 
loops. 

The second input is as ~ g^/(47r), the running gauge 
coupling. It runs with the physical cut-off scale fcphys de- 
rived in AppendixICl as ~ Q^s (fcphys)- present work 
we model with a temperature-dependent coupling that 
runs into a three-dimensional fixed point a*,3dfcphys/2^ for 
low cut-off scales fcphys/^ ^ I. This choice carries some 
uncertainty as the running coupling in Yang-Mills theory 
is not universal beyond two loop order. Here we have 
chosen the Landau gauge couplings aLandau 4d(fc ph vs) 
cut-off scales fcphys/T > 1, see [1, d, [H, il [13] . The 
corresponding three-dimensional fixed point a^.sd = 1.12 
is obtained from [23l | . A specific choice for such a running 
coupling is given in Fig. O 



a. 

til 

a 






MeV 




1 50 MeV 


T = 


300 MeV 


T = 


600 MeV 



0.01 



0.01 



0.1 1 

p [GeV] 



10 



FIG. 2: as for temperatures T = 0, 150, 300, 600 MeV 

The normalisation of the momentum scale has been done 
by the comparison of continuum Landau gauge propaga- 
tors to their lattice analogues. Fixing the lattice string 
tension to ^Ja = 440 MeV, we are led to the above 
momentum scales. For a comparison with the Landau 
gauge results obtained in [l^ we have also computed 
the temperature-dependence of the Polyakov loop by us- 
ing aLandau, 4(i for all cut-off scales. Indeed, this over- 
estimates the strength of as , as can be seen from Fig. [51 
However, qualitatively this does not make a difference: 
for infrared scales far below the temperature scale, fc — > 0, 
the flow switches off for fields (p with d'^{Vj_ + AV) > 0, 
that is for the convex part of the potential. This hap- 
pens both for — > const, and for g^{k^ — > 0) ~ fc. In 
other words, the minimum of the potential freezes out 
in this regime. For the non-convex part of the potential, 
d'^{V± + AV) < 0; the flow does not tend to zero but sim- 
ply flattens the potential, thus arranging for convexity of 
the effective potential Ves = Vk=o- The uncertainty in 
gl is taken into account by evaluating the limiting cases. 
Together with the error estimate on the physical cut- 
off scale fcphys in Appendix [C] this leads to an estimate 
for the systematic error of the results presented below. 



This error includes that related to our specific choice of 
the running coupling. For example, a viable alternative 
choice to Fig. [2] is provided by the background field cou- 
pling derived in [3C| which is covered by the above error 
estimate. 



VI. INTEGRATION 

The numerical solution of ([5^ is done on a suitably 
chosen grid or parameterisation of AV^ and its deriva- 
tives. As V, V± and AV^ are periodic, one is tem pted 
to solve (|32)) in a Fourier decomposition, see e.g. [27| . 
However, as can be seen already at the example of the 
perturbative Weiss potential Vw = Vj^,o, (|B6p . this pe- 
riodicity is deceiving. The Weiss potential is polynomial 
of order four in (p ~ ip mod 27r, its periodicity comes 
from the periodic <^(</3), p^ . Consequently the third 
derivative d^Vw jumps a,t (p = 2?™ with n S Z. More- 
over, d^Vwiip^ ^ 0+] = -d^Vwlv ^ 0_] ^ 0. A pe- 
riodic expansion of Vw, e.g. in trigonometric functions 
cannot capture this property at finite order. This does 
not only destabilise the parameterisation, but also fails 
to capture important physics: the flow of the position of 
the minima is proportional to d^V. This follows from 



dtV[ipr, 



dtip, 



0. Expanding this identity leads to 



'min.k 



V"[^] 



(33) 



v=v>„ 



where V' = d^V and V" = d%V . The flow dtV'{p)\ is pro- 
portional to d^V , which e.g. can be seen by taking the 
(^-derivative of ([19]). Hence, as a Fourier-decomposition 



enforces d^V = at any finite order, the minimum does 
not flow in such an approximation, and the theory always 
remains in the deconfined phase. Note also that the re- 
sulting effective potential at /c = for smooth periodic 
potentials and flows vanishes identically as it has to be 
convex. In the present case this is not so, as the poten- 
tial is rather polynomial (in (p) and convexity does not 
enforce a vanishing effective potential. 

In turn, a standard polynomial expansion about the 
minimum Pmin,fc already captures the flow towards the 
confining phase. Here, however, we use a grid evaluation 
of the flow of Ay with (/? e [0 , 27r] while taking special 
care of the boundary conditions at = 0, 27r: we have 
extrapolated the second derivative to 93 = and (p> — l-n. 
It suffices to use a first order extrapolation, and we have 
explicitly checked that the resulting fiow is insensitive to 
the precision of the extrapolation. 

An alternative procedure is an expansion in terms of 
Chebyshev polynomials that also works quite well and is 
also a very fast and efficient way of integrating the flow. 
A comparison between the results obtained on a grid and 
with Chebyshev polynomials shows that both parameter- 
isations agree nicely and the corresponding flows deviate 
from each other only for small values of fc. This is due to 
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an expected failure of the standard Chebyshev-expansion 
for those small k where the position of the minimum is 
almost settled and the potential flattens out in the re- 
gions that are not convex. This is better resolved with a 
grid than with polynomials. On a grid implementation 
we see the potential becoming convex as fc ^ 0. 



VII. RESULTS 

In Fig. [3] we show the full effective potential for tem- 
peratures ranging from T = 500 McV in the deconfined 
phase to T = 250 MeV in the confined phase. The 
expectation value (ip) in the center-broken deconfined 
phase is given by the transition point between decreasing 
part of the potential for small (f and the flat region in 
the middle of the plot. It can also be explicitly computed 
from (j33p . In the center-symmetric confined phase it is 
just given by the minimum at = tt. 




T = 500 MeV 
T = 400 MeV 
T = 350 MeV - 
T = 300 MeV 
T = 250 MeV 



2 II 



FIG. 3: Full effective potential Kff, normalised to at = 

The temperature-dependence of the order parameter 
L[(Ao)] = cos(((^/2)) is shown in Fig. [H and we observe 
a second order phase transition from the confined to the 
deconfined phase at a critical temperature 



SOSlgg MeV, 



0.691:?^ , 



(34) 



with the string tension ^/a = 440 MeV. The correspond- 
ing value on the lattice is T^j \fa = .709, and agrees 
within the errors with our result. The estimate of the 
systematic error in (|34|) is dominated by that of the un- 
certainty of the identification of fcphys, see Appendix [Cl 

We would also like to comment on the difference of 
the temperature-dependence of _L[(Ao)] depicted in Fig. [4] 
and that of the Polyakov loop (i[Ao]). It has been shown 
in section |TT] that in the confined phase they both vanish 
and both are non-zero in the deconfined phase. However, 
the Jensen inequality jS]) entails that the present observ- 
able L[(Ao)] takes bigger values than the Polyakov loop 
(i[Ao]), which is in agreement with lattice results. 



0.8 



0.6 



0.4 



0.2 




zero temperature coupling 
finite temperature coupling 



1.05 1.1 
T/T 



1.15 



1.2 



FIG. 4: Temperature dependence of the Polyakov loop 
L[(Ao>] = cos({v5)/2) in SC/(2) 



The critical physics should not depend on this issue. 
Here we compute the critical exponent ^, a quantity 
well-studied in the 0(1) model which is in the same 
universality class as SU{2) Yang-Mills theory. More- 
over, in Polyakov gauge the effective action T[Ao] after 
integrating-out the spatial gauge field is close to that of 
an 0(l)-model. Studies using the FRG in local poten- 
tial approximation with an optimised cut-off for the 0(1) 
model yield v = 0.65, see [23. The critical exponent is 
related to the screening mass of temporal gauge field by 



m^{T) oc \T-Tc\ 



(35) 



where — V" {(pmin.o)/2. We have computed the 
temperature-dependence of the screening mass in the 
confined phase near the phase transition, and extracted 
the critical exponent v from a linear fit to the data. This 
is shown in Fig. [5] The fit yields the anticipated value of 



(36) 



for the critical exponent ly. The critical exponent j3 agrees 
within the errors with the Ising exponent (3 = 0.33. 

Finally we would like to compare the results obtained 
here with the results of [l^ . There, the effective potential 
V;3fr[^o] was computed from the flow 0, [U of Landau 
gauge propagators [H, [2^, [2^ within a background field 
approach in Landau-DeWitt gauge. In this gauge the 
confining properties of the theory are encoded in the non- 
trivial momentum dependence of the gluon and ghost 
propagators. Indeed, in [l^ the effective potential 14 was 
computed solely from this momentum dependence and 
was not fed back into the flow. In SU (2) Landau gauge 
Yang-Mills this is expected to be a good approximation 
with the exception of temperatures close to the phase 
transition, see [l^. The back-reaction of the effective 
potential is particularly important for the critical physics, 
and the value of the critical temperature [3l|. 

For the comparison we have computed the present flow 
with the zcro-tcmperature running coupling in Fig. [2] for 
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I I I I , , I I I 

-0.02 -0.04 -0.08 

(T-T,)/T, 

FIG. 5; Critical exponent u from = V" {ipmin,o)/'2 



all temperatures. This mimics the approximation used 
in p^ . which implicitly relies on the zero-temperature 
running coupling as- We also remark that the quan- 
tity L[(Ao)] in general is gauge-dependent, and only the 
critical temperature derived from it is not. However, in 
Landau-DeWitt gauge with backgrounds in Polyakov 
gauge temporal fluctuations about this background in- 
clude those in Polyakov gauge. For this reason we might 
expect a rather quantitative agreement for the quantity 
L[(y4.o)] in both approaches. The results for the tem- 
perature dependence of the Polyakov loop are depicted 
in Fig. [SI The coincidence between the two gauges is 




1 1.05 1.1 1.15 1.2 



FIG. 6: Comparison of L[{ylo)] computed in Polyakov gauge 
and in Landau-DeWitt gauge from [l9l |. 

very remarkable, particularly since the mechanisms driv- 
ing confinement are quite different in the different ap- 
proaches, as are the approximations used in both cases. 
This provides further support for the respective results. 
It also sustains the argument concerning the lack of gauge 
dependence made above. The quantitative deviations in 
the vicinity of the phase transition are due to the trunca- 
tion used in [l^, that cannot encode the correct critical 



physics yet, as has been already discussed there. 



VIII. SUMMARY AND OUTLOOK 

In the present work we have put forward a formula- 
tion of QCD in Polyakov gauge. We have argued that 
this gauge is specifically well-adapted for the investiga- 
tion of the confinement-deconfinement phase transition 
as the order parameter, the Polyakov loop expectation 
value (L[Ao]), has a simple representation in terms of 
the temporal gauge field. Moreover, we have shown that 
L[(Ao)] also serves as an order parameter. In summary 
this allows us to access the phase transition within a sim- 
ple approximation to the full effective action. 

The computation was done for the gauge group SU{2), 
where we observe a second order phase transition at a 
critical temperature of Tc = SOS^gg MeV, as well as the 
Ising critical exponents v and (3 to the precision achieved 
within our approximation. The temperature-dependence 
of the order parameter L[(y4.o)] agrees well with a recent 
computation in Landau gauge [19[ . This is very remark- 
able: firstly the latter computation is technically differ- 
ent as in Landau gauge the full momentum-dependence 
of the propagators is needed to cover confinement. Sec- 
ondly the order parameter L[{Aq)] is gauge dependent, 
only the critical temperature is not. 

In the present analysis we used several external inputs 
which we plan to remove in future work. First of all we 
proceed with computing the running coupling within 
Polyakov gauge, that is the momentum-dependence of 
the temporal gauge field. As it is one of the advan- 
tages of the computation in Polyakov gauge that the 
momentum dependence of Green functions is rather 
mild we expect only minor deviations from the compu- 
tations shown here. As argued in the present work, the 
momentum-dependence of the Ao-propagator also gives 
access to the string tension. For a description of spatial 
confinement one has to treat the spatial components 
of the gauge field beyond the present perturbative ap- 
proximation. Moreover, the present analysis is extended 
to SU{3), which is conceptually straightforward but 
technically more challenging. For the matter sector one 
can revert to the plethora of results with the present 
renormalisation group methods, ranging from results in 
effective theories to that in QCD-based approaches, see 

e.g. i EE ail. 
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APPENDIX A: FADDEEV-POPOV 
DETERMINANT 

From the gauge fixing fimctionals (fT2|) and ([T3|) we can 
compute the Faddeev-Popov determinant given by 



Afp[A] = det 



(Al) 



where A" is the gauge transformed gauge field A. For 
infinitesimal gauge transformations it is given by 

A- = - {d^a^ + ^g4[cT^ d"])^'^ . (A2) 

In the following we use the representation w°cr° = 
cj+(T^ + cj^(T+ + uj^a^, and the related derivatives w.r.t. 
lu^,llj'^. The matrix elements related to w-derivatives of 
F+ read 



6F+{A'^) 
SF+{A'^) 

6F+{A'^) 

5uj3 



-Trcr+ {do(J- +iAl[(T-,<j^]) 



= 0, 



-Ti-a+ {doa^ + iA+[<j^,<7-]) . (A3) 



J 



Analogously we get for the ^-derivatives of F 

SF-{A'^) 



Suj+ 
SF-iA'^) 



0, 



= -Ti-a- {doa+ + iAl[a+,a'']) , 
= -Tr a- {doa^ + lA^ [a\a+]) . (A4) 



Suj 
SF-{A'^) 



The w-derivatives of F^ yield long expressions, and we 
only display the parts proportional to doTia^Ao, where 
we have abbreviated additional terms proportional to the 
spatial gauge fields by dots, 



SF'\A'^) 
Slu+ 

SF^{A'^) 
Suj~ 

SF^iA'^) 
5lo^ 



= -idoA-^Tia^la- ,a+] 



(A5) 



Evaluating the traces (jA3p . (|A4[) . (jA5[) we can compute the Faddeev-Popov determinant. Again we only concentrate 
on the terms dependent on Aq, and use the gauge fixing condition Aq = Aq = for eliminating some of the 
off-diagonal elements. 



-det 



do + igAl 
do- igAl 
~Aig J dxodiA^ H Aig J dxodiAf H l/2(5g + / dxodi H ) 



-dct[{do + igAl){do - igAl)^ (^d^ + J dxodf + J dxodxidi + J dxodxidx2di^] 



(A6) 



Using the third gauge fixing condition, OqAq = 0, we can 
write the Faddeev-Popov determinant as 



^FP 



[A] 



-det 



{d', + {gAiy)]dct[{d', + ■■■)] 



We note that the second determinant in (|A7P is indepen- 
dent of the gauge fields and hence can be absorbed in the 
normalisation of the path integral. The first determinant 
is evaluated in frequency space, we get 



n (5^o(^)) n {i27rTnf~igAUx)f 



(A7) 



Multiplying the determinant (|A7|) with a further constant 
normalisation 



we arrive at 



AAdet 



[G..i=nM^(.))^"ff(i-(^)") 



(A9) 

Eq. (|A9[) is just a product representation of the sine- 
function, sin(a;) = x nn=i° ~ (j^'n)'-' ) ' ^'^^ final re- 
sult for the Faddeev-Popov determinant is 



AFp[A]=M'i2Tf 



n 



sm 



2T 



, (AlO) 



^ = n C^^Tn f 



(A8) where Af' is a further normalisation constant. 
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APPENDIX B: INTEGRATING OUT SPATIAL GLUONS 



After integrating out the longitudinal gauge fields the action ^cff j fx reads 



d*xAf 



I 



A',+OiA!) 



(Bl) 



Writing = ip/(gj3) + a^, where Lp is a. constant and 
ao the fluctuating field, this expression is given to second 
order in the fluctuating fields by 

Sym « \ J drd^x {Zoidaof ~2^r^'{doAf)A^+ 
A- i^{d', + d')6,, - d,d,) A-} 



1 



drd'^ 



[{daof-And^-d,dj)A';- 



AtD^.^D^.'Al} 

where we have defined 



(B2) 



(B3) 



In the present work we neglect back-reactions of the 
potential on the transversal gauge fields. Assuming an 
expansion around Af = 0, F^^^ is block-diagonal, like the 
regulators, cf. eq. (|23|) . and we can decompose the flow 
equation (|19p into a sum of two contributions, schemati- 
cally written as 



'ti k 




Tidt 



dtRk 



ln(4t 



Ra) ..(B4) 



The first term on the rhs encodes the quantum fluctua- 
tions of Aq, the second line encodes those of the transver- 
sal spatial components of the gauge field. In the present 
truncation the second line is a total derivative w.r.t. t, 
and does not receive contributions from the first term. 
Therefore we can evaluate the flow of the second contri- 
bution, and use its output V±,k{AQ) as an input for the 
remaining flow. 

The computation is done for the regulators (j23|l . As 
explained below in scction lTVl the cut-off parameters 
fc, and kj_ in Rk for the fluctuations of Aq and Rk,i_ for 
the fluctuations of A± respectively satisfy a non-trivial 
relation kj_ = k±(k) for coinciding physical infrared cut- 
offs fco for ^0 and k± for A± . The computation is similar 
to those done in one loop perturbation theory in SU (2) 
by Weiss [l^] , the only difference being the infrared cut- 



off. Wc infer from the second line in (jB4p that 

k 

V±,k = l^±.A„v + ^Tr Ihs'^I + Ra)] 



Vw + TY^J^,eikl-p^)Hkl 



(B5) 



Dl). 



In (|B5|) we have used that V^_l,Auv-»oo = up to a con- 
stant term, and have added and subtracted the Weiss 
potential Vw [23 |. 

Vwi^) = -(^ - ^)V(6/3") + - 7r)V(127r2/34) , (B6) 

with the dimensionless ip ~ g/3Ao, and ip = ip mod 2tt. 
Alternatively one can simple put Auv = 0, even though 
this seems to be counter-intuitive. We also have used that 
with it follows trllj^ = 2. Performing the Matsubara 
sum and neglecting terms independent of the temporal 
gauge fields, the resulting effective potential is given by 

^i.fc = dw'{ln(l-2cosMe-^'=- (B7) 



+e 



^20kj 



- ln(l - 2cos(v?)e-'^P + 



From (|B7P we deduce that the potential 14^ approaches 
Vw in the limit k ^ and vanishes like e~^^^ cos((/3) for 
A: — !■ oo. From eq. (jB4[) wc can now extract the flow of 
the effective potential, by setting V^^.k — AV^ 4- 
Then we get 



dtm-. = I 



d^p {7jo{k^ - f) + 2P)9ik^ - 



2 J (2^)3 k^+glf32(AVl' + Vl,) ' 

(B8) 

with the input V±^k given in (|B7|) and 770 = 9* IuZq. The 
factor arises from the fact that we parametrise the 
potential in terms of (p rather than in Aq, and g'^ = g^ /Zq 
is nothing but the running coupling at momentum ~ 
^phys- Thus we estimate gl = Anasiir^ — ^phys)- Note 
that gk is an RG-invariant. The momentum integration 
can be performed analytically, and we are led to 



PdkAVk 



(l + 7]o/5)fc' 



dliV^,k + AVk) 



3(27r)2 1 

where 770 is given by 

Vo = -dtlogas , 
as the consistent choice in the given truncation. 



(B9) 



(BIO) 
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APPENDIX C: MATCHING SCALES 

The flow of the temporal component of the gauge field, 
Aq(x), is computed with a three-dimensional regulator, 
see In Polyakov gauge Aq{x) only depends on 

the spatial coordinates, whereas the spatial components 
A±{x) are four-dimensional fields. For cut-off scales far 
lower than the temperature, k/T <C 1, also the spa- 
tial gauge fields are effectively three-dimensional fields 
as only the Matsubara zero mode propagates. Hence in 
this regime we can identify k = k±. For large cut-off 
scales, k/T 3> 1, the Ag-flow decouples from the theory. 
A comparison between the two Sows can only be done 
after the summation of the spatial flow over the Matsub- 
ara frequencies. In the asymptotic regime k/T 3> 1 this 
leads to the relation 

1 °° 1 1 

fc~ ^ w2 ^ 2fcl ' ^^^^ 

n— — oo 

The crossover between these asymptotic regimes happens 
at about k/T = 1. This crossover is implemented with 
the help of an appropriately chosen interpolating function 

n— — oo 

A natural choice for f{k/T) is depicted in Fig. [71 and 
has been used in the computation. A more sophisti- 
cated adjustment of the relative scales can be performed 
within a comparison of the flow of momentum-dependent 
observables such as the wave function rcnormalisation 
Zq. The peak of these flows in momentum space is 
directly related to the cut-off scale. Indeed, the func- 
tion / carries the physical information of the peak of 
the flow at some momentum scale. Scanning the set of 
/ gives some further access to the uncertainty in such 
a procedure. The effective cut-off scales fcphys(fco) and 




0.5 L 



5 10 15 20 

A 

k 

FIG. 7: k_i_/k as function of k. 



fc±,phys(fcj_) in the flows of the temporal gluons and of 
spatial gluons respectively do not match in general. If 
solving the flow within a local truncation as chosen in 
the present work we have to identify the two effective 
cut-ofl: scales, fcphys(fco) = fc_L_phys(fcj_) = A:phys, leading 
to a non-trivial relation fco = ko{kj_). Moreover, the ef- 
fective cut-off scale has to be used in the running coupling 
as = as{f = fc^j^yj. 




0.6 I • \ • \ • \ • 1 

5 10 15 20 

k 

FIG. 8: fcphys(fc) from the comparison of flows with three- 
dimensional regulators and four-dimensional regulators. 



It is left to determine the physical cut-off scale /cphys 
from cither the flow of the spatial gauge fields as 
fci,phys(fcj_) or from the temporal flow fco,phys(fco)- We 
first discuss the spatial flow. For an optimised regulator 
depending on all momentum directions, , we have the 
relation fcphys = fcj.. Hence the relation A:^.phys(^±) can 
be computed if comparing the flows for a specific observ- 
able with three-dimensional regulator Ropt,k±ip^), (EH), 
with fiows with four-dimensional regulator -Ropt, fcphys (p^)- 
Here, as a model example, we choose the effective poten- 
tial of a (/(^-theory. This leads to the relation fcphys (fcj.) 
displayed in Fig [51 We remark that the relation in 
Fig. [51 depends on the dimension d of the theory, and 
flatten to fcphys(fc) = k for d — > oo. In other words, 
limfc^oo fcphys (fc)/fc is proportional to d/{d— 1). More- 
over, for momentum-dependent observables the crossover 
rather resembles the relation fc_L(fco) as it is more sensi- 
tive to the propagator than to the momentum integral of 
the propagator. Indeed, for the three-dimensional field 
Aq(x) the cut-off scale fco is another natural choice for 
the physical cut-off scale, fco,phys(fco) = fco, even though 
it underestimates the importance of the spatial flow for 
the correlations of the temporal gauge field. In summary, 
we take the above two extremal choices fcphys = fc±,phys 
depicted in Fig.[51and fcphys = fco,phys(fco) = fco as a broad 
estimate of the systematic error in the present computa- 
tion. 
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